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Numerical Simulations of (f = 3 SU (2) LGT in the Dual Formulation 
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We have developed the techniques necessary for a numerical simulation of d — 3 517(2) Lattice Gauge theories 
in the dual formulation as originally developed by Anishetty et al These include updating techniques that 
preserve the constrained configuration space, efficient evaluation of 6-j symbols and a certain problem associated 
with the positive indefiniteness of the weight factors. 



1. Motivation 

The outstanding problems in QCD are un- 
derstanding the mechanisms of confinement, chi- 
ral symmetry breaking and the spectrum of the 
theory including the vacuum structure. Lattice 
gauge theories have gone a long way in shedding 
light on some of these important issues. Never- 
theless, we are still not in a position to, say, dis- 
criminate between different mechanisms for con- 
finement like the dual superconductor mecha- 
nism, the Z/v-fluxon mechanisms etc. In particu- 
lar a detailed picture of QCD in terms of the dual 
variables is still at a preliminary level. Therefore 
it is desirable to probe alternative formulations 
of LGT that place greater emphasis on the dual 
variables. 

The dual variables simplify the description of 
some observables while complicating some oth- 
ers. In the formulation to be discussed here, we 
shall see that while the t'Hooft disordered vari- 
ables become diagonal, the Wilson loop operators 
become quite complicated. The continuum limit 
of LGT's is in the extreme weak coupling regime 
and here many of the gauge invariant observables 
like plaquettes approach unity so in extracting the 
scaling part of various correlation lengths large 
"perturbative" corrections need to be subtracted. 
But in the dual formulation there are many gauge 
invariant observables that actually vanish in the 
weak coupling regime so it can be hoped that 
correlation functions of these observables are lot 
"cleaner". Other important feature of the dual 
formulation is that it provides manifestly gauge 



invariant characterisation of fluctuations. 

As shown by Q there is an intriguing connec- 
tion between the dual formulation of Yang-Mills 
theories in d = 3 and the Regge-Ponzano formu- 
lation of gravity in d = 3. A count of the relevant 
degrees of freedom indicates that this should hap- 
pen in d = 4 also. 

2. The Dual Formulation 

In the conventional formulation, group ele- 
ments are assigned to the links of the lattice and 
gauge invariant action is composed of the traces 
of the plaquette variables with the partition sum 
given by 
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In order to go over to the dual formulation one 
introduces the group character expansion 

exptftrfPi) = £(2j< + l)C jt ([3) Xjl (Pi) (2) 



Now either group integrations over products of ro- 
tation matrices or graphical rules j^j can be used 
to carry out the link integrations. In d = 3, every 
link is shared by 4 plaquettes and the link integral 
is nontrivial if the representations attached to the 
four plaquettes are such that ji®j2®j3®ji has a 
singlet in it. By introducing an auxiliary j§, this 
condition can be stated as (ji, j 2 , js), (j3, J4, js) 
should satisfy the triangle inequalities. These are 
the analog of Gauss's law constraints of the con- 
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ventional formulation. For lack of space, I'll just 
quote the final result |Q 

z d =^m+mMi[{2 % 2} (3) 

This partition sum is defined over the dual lat- 
tice where {j a } live over the links and {jb} over 
the diagonals to the plaquettes. The convention 
for the diagonals is that they connect the vertices 
of the odd sublattice.We have collectively desig- 
nated {j a },{jb} by {j}. Each cube of the dual 
lattice is seen to be spanned by 5 tetrahedra of 
which one is spanned entirely by {jb} while four 
are spanned by three {j a } and three {jb}- Each 
tetrahedron carries a weight factor which is the 

SU{2) 6-j symol / ^ ^ ^j. Periodic b.c for 

the original lattice is crucial for this construction. 

The disordered operators Di satisfy 
X n { p i) {-^XhiPi)- The disordered line 
is composed of the product of disordered opera- 
tors for each plaquette pierced by the disordered 
line. 

2.1. Connection With Regge-Ponzano 
Gravity 

The partition function of eqn (3) is the same 
as the partition function of the d = 3 Regge- 
Ponzano gravity except for the Cj a factors. These 
factors break the invariance under the Alexander 
moves |lj thereby generating an additional de- 
gree of freedom per generator of the gauge group. 
SU(2) gauge theory in D dimensions has 3(1? — 2) 
d.o.f. The number of diffeomorphisms being D, 
the d.o.f left is 2D — 6. This correctly matches 
the d.o.f of gravity in D = 3, 4 dimensions. 

3. Numerical Simulations 

Apart from the new feature of simulating in 
the dual variables, there are many practical ad- 
vantages also. The basic variables are all inte- 
gers. Further, since Cj{[3) e ~ 3 =0 , only 
jmax — VP is needed making short integers to 
suffice. 



3.1. Updating 

An allowed configuration is where various tri- 
angle inequalities are satisfied and any updating 
of variables has to maintain this. Updating {j} 
independently and rejecting configurations that 
are not allowed is wasteful. First we describe 
what we call a local updating algorithm. Each 
a-link is common to 4 triangles and each b-link is 
common to 6 triangles. To update a link, one first 
forms j max = j\ 4- j\ for each triangle i contain- 
ing the link and where j± t 2 are the reps on the 
other two links of the triangle. Likewise j min = 
\jl ~ ill are a lso calculated and subsequently 
jmax = min{f max } and j min = max{j l mln }. Now 
the allowed values for the link to be updated are 
in the range (j m i n ,jmax) in steps of unity. A par- 
ticular value can be chosen by either the heatbath 
or the Metropolis algorithms. 

A major problem with this updating is that 
it does not change a link from integer val- 
ues (bosonic) to half- integer values (fermionic) and 
vice versa. Therefore it is not ergodic. If we 
think of the integer and half-integer values as a 
/^-grading, then these Zi d.o.f are not updated 
by the local method. 

Changing only one link from fermionic to 
bosonic or vice versa is not possible. In any 
given triangle at least two such changes must be 
made. Continuing this way one is generically led 
to a proliferation of changes! We must seek the 
smallest volume to which such changes can be re- 
stricted. If the Z2 flippings are carried out at an 
odd site, the minimum volume is 8 cubes. These 
flippings called Quasi-local updates amount to 
keeping all the links on the outer surface fixed 
and flipping all the interrior links. If on the other 
hand the flippings are done at an even site, the 
minimum volume is a cluster of 8 tetrahedra. 

One still faces the question of how these flip- 
pings are to be done satisfying all the trian- 
gle inequalities. That brings us to the use of 
Kagome(K) variables ||. The basic idea is that 
the triangle inequalities can be traded for equal- 
ities. If (31,32,33) are the links of a triangle, 
introducing the variables j\ = ?i2 + ^3 , 31 = 
n 3 + n i , 33 — ni+ri2, it can be shown that choos- 
ing m freely among positive integers, one can gen- 
erate all the allowed ji. But for our case, if a link 



3 



is shared between two triangles, then the sum of 
K-variables for this link from one triangle must 
equal the sum of the corresponding K-variables 
for the other triangle. Thus K-variables are not 
all independent. 

Now the solution to the quasi-local update 
problem in terms of the K-variables is as fol- 
lows: i) determine the K-variables for the links 
on the outer surface and keep them hxed; ii) de- 
termine the the set of independent K-variables for 
the interrior; iii) change the interrior variables 
freely subject to realising the fermion-boson flip- 
pings for the interrior links. A simplified quasi- 
local update consists of simply shifting all {j} by 
{J ±1/2}- 

The ^-fluctuations can be controlled by 
switching on/off the quasi- local updates at any 
site one wishes. Thermalisation can be achieved 
first with the local updates and on the resulting 
ensemble one can perform the quasi-local updates 
till the ^-configurations thcrmalisc. In this man- 
ner the precise influence of the ^-fluctuations on 
various dynamical mechanisms can be studied. 

Finally it is worth mentioning as to how best 
the 6-j evaluations should be done. Calculating 
them for each configuration is unnecessarily time 
consuming. For local updates it is possible to use 
a recursive algorithm. But it is best to make a 
table of only the non- vanishing 6-j's into an one- 
dimensional array with an unique labelling. 

3.2. Positivity of Weights? 

The product of 6-j symbols is generally not of 
a particular sign. Nevertheless, based on my ex- 
periences with the d = 4 hypercube problem |2j 
it was initially hoped that a configuration of pos- 
itive weight remains so under all updates. Over 
a large number of sweeps this was indeed so, but 
eventually negative weight configurations started 
appearing. Barring a phase error in the analytical 
calculations, it appeared as if the dual formula- 
tion had seen the end of the line at least as far as 
numerical simulations were concerned. 

But a resolution of this problem has recently 
occurred to me. Let us say the weight factors 
Pi,P2tP3) ■■■■iPki ■■■ are such that some of them, 
say, PitPmiPn, ••• are negative. Now the trick 
is to generate configurations with the probabil- 



ities {pi — -^r-} where P — Yli \pi\ but keep 
track of the configurations with negative weight 
through €i which is +1 when pi > and —1 when 
Pi < O.The expectation value of an observable 
with value Oi in the configuration i is now given 
by 



(0) P = 



(Oe)p 



Since the partition sum (l) p 
procedure is well defined. 

4. Open Issues 



(4) 

<e)p ^ 0, this 



Among the important open problems are the 
generalisation of our algorithms to d = 4 us- 
ing the results of It would also be interest- 
ing to generalise our results for the case of the 
models considered by Ooguri[^| and their non- 
topological versions. Inclusion of matter partic- 
ularly fermions is another outstanding problem. 
Since we are eventually interested only in the 
large [3 regime a simpler version using only the 
asymptotic forms of the 6-j symbols may be pos- 
sible. Generalisation to SU(3) is another open 
direction. Finally mention must be made of what 
I call exceptional configurations for which some 
of the 6-j's vanish even though all the triangle 
inequalities are satisfied. There are fascinating 
group theoretical explanations for these. Their 
contribution to the partition sum is zero. One 
may speculate as to the physical significance of 
these configurations for gauge theories. 
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